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Abstract 



Using the new variational approach proposed recently for a systematic improve- 
ment of the locally harmonic Feynman-Kleinert approximation to path integrals 
we calculate the partition function of the anharmonic oscillator for all tempera- 
tures and coupling strength with high accuracy. 

1) Some time ago, Feynman and Kleinert [1] extended the variational approach to eu- 
clidean path integrals developed earlier by Feynman in his textbook on statistical me- 
chanics [2], improving greatly the accuracy at low temperatures. A similar extension 
was found independently by Giacetti and Tognetti [3] who applied this method to sev- 
eral statistical systems [4]. The approach is particularly successful in systems in which 
the quantum effects do not produce essentially new phenomena, their main result being 
a modification of the quasi-harmonic properties of a system. This happens in quantum 
crystals [5] where experimental data can now be explained very well. The approach 
yields a good approximation to the effective classical potential Kff,ci at all temperatures 
and serves to calculate the free energy of the system as well as particle distributions at 
all coupling strength, including the strong coupling limit. 

The purpose of this note is to present results of a recently proposed systematic 
improvement [6] of this method for the anharmonic oscillator carried out to third order 
in the coupling strength. The results turn out to be in excellent agreement with a 
precise numerical values of the free energy [7] at all temperatures and coupling strengths 
available in the literature. 

2) Let us briefiy review the Feynman-Kleinert approach. The aim is to write the parti- 
tion function as a classical configuration space integral over a Boltzmann factor involving 
an effective classical potential VeH,d{xo): 
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with /3 = l/k^T . The variable Xq is equal to the time-averaged position x = J^^^ dTx{T) 
of the fluctuating path. The effective potential Kff.cil^^o) is defined by the path integral 

exp {-/3Kff,ci} = = jjjVx 8{x - Xo) exp {-Ajh} , (2) 
with the euclidean action 

A= I dT{ -Mi^ + V{x] 

where 



^2 



S{x — Xq) = ■\J 27^%^ (3 1 M6{x — Xq) 



restricts the path average x to the value Xq. The paths are periodic in fifi. We shall 
refer to as a restricted partition function . The usefulness of separating out x derives 
from the fact that at any temperature the fiuctuations rarely carry x{t) far away from 
X. This allows approximating the deviations from x in a quasi-harmonic way. One 
reexpresses the local partition function as an expectation value within a harmonic trial 
system whose action is, for each Xq, 

M f^^ 

= V / + - ^o)'} (3) 
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and chooses VP{xo) optimally. With the help of , the defining expression (2) can be 
rewritten as 

exp{-/3Kfr,ciM} = {Z^^^V {ew{-{A- A^^)lh}Y^ (4) 

for any 0(a;o). The subscript stands for the trial system in which the expectation is 
calculated and the superscript Xq indicates the restricted value of x. 

The expectation value on the right-hand side of (4) cannot be calculated exactly. A 
lowest-order approximation is found with the help of the Jensen-Peierls inequality 

(exp {-{A - Al^)in}Y^ > exp {- - Air^ Ih] . 

The harmonic expectation value in the exponent of the right-hand side yields the ap- 
proximation Wi for the effective classical potential 

V^sA^o) ~ W^ixo) = V^^ixo) + V,2{xo) - ^Mn'{xo)a\xo). (5) 

The first term is the logarithm of the restricted partition function of the trial system 
— In {(/i/30(a;o)/2)/ sinh(/i/30(a;o)/2)} the last two terms are the restricted expecta- 
tions 

{V{x)rn° = VAxo) 
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and 

\m^\x,) {{x - x,fY^ = \M^\x,)a^{x,). 
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In the last equation we have abbreviated the restricted square deviation ((x — XoYY° 



The first term is the usual thermal expectation < x'^{t) >; the second term subtracts 
from this the square deviations of the temporal average x from Xq. The restricted 
expectation of the potential is obtained by a simple gaussian smearing process of a 
width a!^{xo) 

-oo V/!vra^ 

The best approximation is obtained by minimizing the function W\{xo) with respect 
to 0(a;o) which yields the condition 

The resulting W\{xo) is always slightly larger than the exact effective classical po- 
tential Kff,ci(a;o)- 

Eqs. (6) and (8) are solved numerically. The resulting approximation to the partition 
function 

dxQ 



Zi = / e 

leads to the free energy Fi = —k^TXnZi which describes the true free energy F = 
— k^TXnZ of the system quite well at all temperatures. Ref. [8] explains why this is so. 

For the anharmonic oscillator potential 

V{x) = —x + -X , 
Eq. (8) leads to the following equation for the optimal 0: 

2 2 , o 1 f .1 



0^ = ?,gxi + 3^/-^ I coth - 1 j + 1 . (9) 

The resulting free energies are listed in Tab. I for comparison with the improved results 
to be derived in this note. 

3) To go beyond the above approximation we split the action into a free part Aq and 
an interacting part ^j^^. The first contains all classical terms depending on Xq plus the 
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trial action A'^ of Eq. (3), the second term all the rest. The partition function can then 
be rewritten as 

Z = I , = r , e^(-")/^ (e-A?.A\^" . (10) 



In the previous section, the Jensen-Peierls inequality was used to approximate the ex- 
pectation value on the right-hand side of (10). Now we evaluate this expectation per- 
turbatively, expanding the exponential function in a Taylor series 

yxo _ -(5V{xo) yxo j 1 _ ^ / A^o \^o , ^ If axq \2\^o _ ^ // axq \3\^o , 

and going over to the connected parts by a reexpansion into cumulants 

Z" = exp{ - /SI'S" - I (-^S,)? + ((-4g,)'>^; 

The truncated exponents 
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•••+^((An"t)T; (12) 

are optimized in ft{xo) and yield the successive approximations to the effective classical 
potential VeH,d{xo). The lowest approximation Wi{xo) obviously coincides with the 
Feynman-Kleinert approximation. 

We study the behavior of Wn{xo) up to = 3. The condition of an optimal ft{xo) 
cannot always be met by an extremum of Wn{xo)^ since for higher > 1 the condition 
dWN{xo) / dft{xo) = does not always have real solutions il[xo). This happens for 
example at = 2. In such cases, an optimal result is obtained by enforcing a minimal 
dependence on via the disappearance of the second derivative of Wn{xo) with respect 
to il[xo). 

To calculate the restricted expectation values in equation (12) we consider polynomial- 
like potentials in view of the intended quantum field theoretic applications. The inter- 
action A^°^ is expanded around Xq 

An"t = f'dr {|(x - XoY + |(x - XoY + |(x - xor • • •} , (13) 
defining the coupling constants 

g, = V^'\xo)-nH,,2. (14) 
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Then (12) is obtained by calculating all connected vacuum diagrams formed with vertices 
gn/n\h and lines which stand for the correlation function 

^ r;^^f^cosh[Ml(l_lbrd)] ] 

in which the first term is the usual Green function {x[t)x[t')) of the harmonic oscillator 
while the second term subtracts from this the zero-frequency part (xq) = l//30^. 

The number of Feynman integrals to be calculated is reduced in three ways. First, 
only one-particle irreducible vacuum graphs contribute (i.e., there are no tadpole dia- 
grams) due to the absence of the zero-frequency modes in 6x[t) = x[t) — Xq. Second, 
Feynman integrals of subdiagrams which touch the rest in one vertex can be factored 
out. Third, diagrams with vertices g2 (mass insertions) can be obtained from those 
without these vertices by replacing the frequency by = ^fl^ -\- §2 and expanding 
the result in powers of §2- 

4) We now turn to the application to an anharmonic oscillator with the potential 

V{x) = —x +-X . 

The interaction is 

An"t = / dT{^{8x)\r) + |(^x)3(r) + |(^x)^(r)} 

(we use natural units with fi = M = 1), with the coupling constants 

g2 = {io' -n') + 3gxl 

93 = Qgxo (16) 

94 = ^9- 

The quadratic coupling g2{Sxy[T) /2 is now included into the harmonic trial potential 
term by changing it to 0^(a; — XqY /2 with the modified frequency introduced above. 
The third-order approximation 

has the diagrammatic expansion 

^3(0) = 




+ ooo 

72 / 
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■ \ 2592 I ) 

648 648 3456 1728 ^"28 

The number under each diagram count in how many Wick contractions it appears; this 
number multiphes the pref actors. 

Only five integrals corresponding to the diagrams 

O ; ^ ; © ; © ; @ 

need to be evaluated explicitly. The remaining diagrams follow by differentiation with 
respect to 0^ and factorization into subdiagrams touching each other at a point. 

In quantum field theory one is used to calculate Feynman diagrams in momentum 
space. At a finite temperature this requires performing multiple sums over Matsubara 
frequencies. To higher order, they are hard to do in closed form. In the present case 
of a D = I -dimensional quantum field theory it is more convenient to evaluate the 
diagrams in r-space as integrals over products of the Green function (15). For example, 
the diagrams 

o, O 

require the integrals (before replacing — )• 0). 

Jo 




respectively. The factors hjS are due to an overall r-integral, the powers of I/O arise from 
the remaining ones. The parameters have the dimension of a length, the subscript 
indicating the number of vertices in the diagram (defining af = a^). In terms of these 
two Feynman integrals, the integral associated with the last diagram in (18) is given by 
the product 
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The symbol = indicates that the right-hand side is only the Feynman integral of the 
diagram without the factors gn/n\h associated with the vertices. 

The explicit expressions for the different Feynman integrals are, with the abbrevia- 
tion X = fifiVl, 



O 



.10 



.12 



Tr log ^(2) = _2j3V^ 



n 



-2 log 



sinh(a;/2) 



( ^ ) «2 



I —24 — 4 -|- 24 cosh x -\- x'^ cosh x — 9 x sinh x 



I ^ ) «2 



1 



1 



nj 768x3 ginh^ 

2 

-288 cosh 2x - 32 x^ cosh 2a; 
-|-f 44 X sinh 2a; -|- 3 x"^ sinh 2a; 



{-864 + f8 + ff52 cosh a; + 32 a;2cosh, 
288 X sinh a; + 24 a;"^ sinh x 



M ( i 



.10 



f 



f 



{ -3456 - 4f 4 x"^ -6 x^ + 4608 cosh x 



n 

n) 2304x^~;^T^ 
2 

+496 a;^ cosh x — f f 52 cosh 2a; — 82 a;^ cosh 2a; — f 008 x sinh x 
— f 6 a;"^ sinh x + 504 x sinh 2a; + 5 a;"^ sinh 2a; } 



.12 







49f 52a;4 gj^^j^e 



{-f 07520 - 7360 a;^ + 624 a;^ + 96 



+ f 6f 280 cosh X + f 2000 a;^ cosh a; — 777 x"^ cosh a; + 24 a;^ cosh x 
— 645f 2 cosh 2a; — 5952 x^ cosh 2a; + f 44 a;^ cosh 2a; + f 0752 cosh 3a; 
+ f 3f 2 a;^ cosh 3a; + 9 a;"* cosh 3a; — 28800 x sinh x + f f 20 x^ sinh x 
+324 x^ sinh a;23040 x sinh 2a; — 320 x^ sinh 2a; — 5760 x sinh 3a; 
-f60 a;3sinh3a;| ff8l 



To obtain the diagrams involving g2 we replace, as explained above, by and expand 
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in powers of (/2/2/i: 



00 = c» - 00 +i cx) +7 00 



24 



72 



144 



108 



1 

2 




24 



24 



576 



1 



72 



, -i ODO -i 

Z D D 

72 864 



COO 

864 



Note that the dot being associated with the vertex (72 /2 is equal to half a mass insertion. 
The numbers underneath the diagrams are their multiphcities; they act as factors. 

The expansion of the trace log follows a different pattern since the one-loop diagram 
requires a reexpansion of the logarithm pictured by it: 



0'= o -20 + <o 4 



As an example, consider the Feynman diagram 

<0 =2(1)^,44 

2 

It is obtained from the second-order Taylor expansion term of the trace log as follows: 



\2h) ' Vt ^ 2! (902)2 



-2/3^5" 



(0^ _ 



2\2 



The differentiation yields 



hV I I 



0/ gJj^]-^2 



|4 -|- — 4 cosh X -\- X sinh x 



corresponding to the Feynman integral 

"'1^ ^ L L 



G'^(t2 — Ti) dTidT2. 
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Similarly we find the Feynman integrals 



( ^ ) «3 




fl\ 1 if X 3 ^ ^ 3x 

— ^ < —3 X cosh \- 2 X cosh \- 3 x cosh — 

J 64x ginh^ - I 2 2 2 



-48 smh \- Q X smh 16 smh — 

2 2 2 



(19) 



( ^ ) «3 



1 



,-y* ,-y* -< ,-y* 

™ , 45 x cosh 6 x cosh 45 x cosh — 

^1 2 2 2 



(20) 



-432 smh 54 x smh — \- 144 smh \- A x smh — 

2 2 2 2 
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/i/? ( - ) Og, 



h\ 1 1 ( X 

— ^ { 672 X cosh — 

y 4096x3 gij^j^5 ^ I 2 



8 cosh - + 24 cosh - 



J rfi -\ rfi l-\ rfi l-\ rv* 

-1008 x cosh h 3 x cosh h 336 x cosh h 5 x cosh — 

2 2 2 2 

3 

-7680 sinh - - 352 x^ sinh - + 72 sinh - + 3840 sinh — 
2 2 2 2 



3^ 3^ 

+224 x^ sinh h 12 sinh — 

2 2 



768 sinh 



5a; 



64 sinh 



5a; 



(21) 



obtain the graphical expansion for W-ii^x^) (the multiplicities of the diagrams have been 
canceled by the n!-denominators of the couplings gnjhn\] now the vertices represent only 

mA^.) = v{xo) - ^0+(^0+^00) 

^c^ +^oo +^ooo 
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+3(iO0O+i©+i00o) 



+ (^0000 +\ 



1 

+ 8 



+ 



corresponding to the analytic expression (in the same order) 
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2mn 
1 



2 222 

92 4 I 9294 4 2 I 5*3 6 I 5*4 8 I 5*4 4 4 

— a, H a^a -\ a, -\ a, -\ a^a 

2^ 2 2 6^ 24^ 8^ 



,5253 8 



2 2 
9294, 4x2 I 5254 6 2 



+ 



+3 



5354 



8 2 I 



5354 



,10 



2 22 

, Q / 5254 / 4\2 2 I 5254 10 I 5254 6 4 



+ 



354 / 4 



16 



2\2 



6 



I 54 10 2 I 54 12 I 54 6 4 

H act; a H a-^ H a-^a 

4 ^ 8^8^ 



f22l 



The best 0(a;o) is found numerically by searching for the roots of the first derivate 
of W3 with respect to 0. There are many solutions and we must pick the right one, 
which we take to be the one closest to the unique solution found at the lowest order 
(where the Jensen-Peierls inequality ensured its existence). 

Fig. 1 shows the approximations VKi,2,3 over the variational parameter for = 4 
and /3 = 1 at Xo = 0. The new extremum lies indeed near the old one. Note that W2 has 
no extremum as remarked before, but the curve has a very small slope and the point of 
smallest 0-dependence is easily identified. In fact, the 0-dependence decreases rapidly 
with increasing order refiecting the fact that the exact Woo is completely independent 
of the choice of 0. 

The corresponding results for the free energy 

1 



F 



■\nZ 



are listed in Tab. 1 and compared with the precise numbers Fgx of Ref. [7]. 

To third order, the new approximation gives energies which are better than the 
lowest order results by a factor of 30 to 50. They differ from the exact results only in 
the forth digit. 
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Figure 1: Effective potential V^fj versus for = 4 and /3 = 1 at = 0. The curves 
1,2,3 show the approximations VKi,2,3 respectively. While Wi and W3 have an extremum 
at 3, 7, the second-order approximation W2 has only extrema in the complex plane. 
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r 


F^ 

' i 




F 

' ex 


0.002 


2.0 


0.427937 


0.427937 


0.427741 


0.4 


1.0 


226084 


226075 


226074 




5.0 


0.559155 


0.558678 


0.558675 


2.0 


1.0 


0.492685 


0.492578 


0.492579 




5.0 


0.699431 


0.696180 


0.696118 




10.0 


0.700934 


0.696285 


0.696176 


4.0 


1.0 


0.657396 


0.6571051 


0.6571049 




5.0 


0.809835 


0.803911 


0.803758 


20 


1.0 


1.18102 


1.17864 


1.17863 




5.0 


1.24158 


1.22516 


1.22459 




10.0 


1.24353 


1.22515 


1.22459 


200 


5.0 


2.54587 


2.50117 


2.49971 


2000 


0.1 


2.6997 


2.69834 


2.69834 




1.0 


5.40827 


5.32319 


5.31989 




10.0 


5.4525 


5.3225 


5.3199 


80000 


0.1 


18.1517 


18.0470 


18.0451 




3.0 


18.501 


18.146 


18.137 



Table 1: Free energy of the anharmonic oscillators with the potential V{x) = 1x^/2 + 
gx'^/A for different g and /3 = 1/kT. 
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Figure 2: The approximations VKi,2,3 for the effective classical potential V^fj versus /3 for 
g = 4: and Xq = with optimized . The curves 1,2,3 shows the approximations VKi,2,3 
respectively. The dotted line is the asymptote of W3. 

In the high-temperature limit, all Wn tend to the classical results (which, in turn, 
becomes exact). It is therefore not astonishing that for small /3, the approximations W3 
and Wi are practically indistinguishable. Fig. 2 shows the dependence of VKi,2,3 on /3. 
For increasing /3, the difference between the first approximation and the other curves 
increases. 

The worst possible case is T = 0, where the diagrams become most simple and have 
been given before [6]. 

1 Conclusion 

We have presented a systematic improvement of the Feynman-Kleinert approximation 
up to third-order. For the anharmonic oscillator the free energy was obtained with an 
error less than 0.04 % for all coupling constants g and all temperatures. In comparison 
with the previously available lowest-order approximation, the accuracy is increased by 
about a factor 30. 

The virtue of the improved variational approach lies in the ability to yield very 
precise results for all coupling constants and temperatures. This is in contrast to other 
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~^ 10 20 30 40 50 ^ 

Figure 3: The effective dimensionless coupling constant g'h/Q,^ = gh/Q,^[g) of the new 
perturbation series versus gh/Q,^ in the first and third approximation Wi and W3, re- 
spectivly (axes are labeled setting h = = 1). 

methods which are applicable only in selected limiting temperature or coupling constant 
regimes. 

The improvement with respect to the ordinary perturbation expansion is due to the 
fact that the new perturbative expansion has an effective coupling constant g'h/Q,^ = 
gh/Q,^[g)^ which remains small even for — )• 00 (see Figure 3). 

Certainly, the accuracy can be increased by calculating higher-order terms. The 
complexity, however, increases rapidly with increasing orders. Thus higher-order calcu- 
lations will require much more work. 

In this note, we have presented only the free energy of the anharmonic oscillator. 
Other quantities like particle distributions can be obtained quite similarly. This will be 
done in future work. 
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